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Abstract 



Theories beyond the Standard Model must necessarily respect its gauge symmetry This 
implies strict constraints on the possible models of non-standard neutrino interactions, which 
we analyze. The focus is set on the effective low-energy dimension six and eight operators 
involving four leptons, decomposing them according to all possible tree-level mediators, 
as a guide for model building. The new couplings are required to have sizeable strength, 
while processes involving four charged leptons are required to be suppressed. For non- 
standard interactions in matter, only diagonal tau-neutrino interactions can escape these 
requirements and can be allowed to result from dimension six operators. Large non-standard 
neutrino interactions from dimension eight operators alone are phenomenologically allowed 
in all flavour channels and shown to require at least two new mediator particles. The new 
couplings must obey general cancellation conditions both at the dimension six and eight 
levels, which result from expressing the operators obtained from the mediator analysis in 
terms of a complete basis of operators. We illustrate with one example how to apply this 
information to model building. 
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1 Introduction 



The experimental observation of neutrino masses and mixings is the first evidence for physics 
beyond the Standard Model (SM) - maybe together with the indication for dark matter - 
and points to the existence of a new, yet unknown, physics scale. The tiny masses of 
the neutrinos, which are orders of magnitude lighter than those of other fermions, suggest 
a large new physics scale leading to very suppressed effects. Since neutrinos have only 
weak interactions with the SM particles, they may even constitute an excellent window into 
the new physics underlying the "dark sectors" of the universe, i.e., dark matter and dark 
energy. Therefore, new physics may very well appear next in the form of exotic couplings 
involving neutrinos, which are often called non-standard neutrino interactions (NSI) [1-5]. 
This possibility is being actively explored and will be the subject of the present work. In 
particular, we will, in a model-independent way, discuss the connection between NSI and 
the possible tree-level mediators of new physics inducing them. If NSI are detected, this 
study will serve as a guide for the model builder. 

Note that the very tiny neutrino mass differences have only been detectable because the 
masses affect the neutrino propagation by inducing small phase shifts, which can be compen- 
sated by the very long distances travelled in neutrino oscillation experiments. In contrast, 
non-standard couplings are short distance (local) effects, which usually do not benefit from 
such an enhancement - unless they affect the propagation in matter, one of the several 
possibilities explored below. Notice, though, that neutrino oscillation experiments may well 
turn out to be the best arena to detect NSI, as they can affect oscillation amplitudes linearly, 
instead of quadratically as in most charged lepton flavour violation processes. 

On general grounds, whatever the nature of the new putative couplings is, observable effects 
will only be expected for a new physics scale A near the present experimental limits, i.e., 
above the electroweak symmetry breaking (EWSB) scale. An example of NSI is given by 
the dimension six (d = 6) operator in 



In this expression, spinor indices are omitted, Greek letters denote flavour indices, Pl = 
(1 — 7s)/2 is the left-handed chiral projection operator, and v and t label the SM neutrino 
and charged lepton fields, respectively. 

The operator in Eq. (1) is not a singlet of the SM gauge group, SU (3) x SU (2) x U(l), while 
the high-energy theory has to contain and encompass the SM gauge group, though. For 
instance, the coupling in Eq. (1) could result from the following gauge invariant operator: 



where L denotes the leptonic SU (2) doublets. Eq. (2) illustrates the consequences of elec- 
troweak gauge invariance: The coupling in Eq. (1) comes necessarily accompanied by other 
charged lepton transitions for which stringent limits may exist. As an example, for (3 = \x 
and a = 7 = 5 = e, Eq. (2) would imply \x — » 3e transitions with the same strength than 
the interaction in Eq. (1). 
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In order to discuss such constraints coming from the gauge invariant framework, it is con- 
venient to rephrase the problem in terms of a generic low-energy effective theory. Effective 
theories allow rather model-independent analyses based on the fundamental symmetries, 
while only the coefficients of the effective operators are model-dependent. The impact of 
the heavy fields present in the high-energy theory can be parametrized, without loss of gen- 
erality, by the addition to the Lagrangian of a tower of non-renormalizable operators O d of 
dimension d > 4, made out of the SM fields and invariant under the SM gauge group. The 
operator coefficients are weighted by inverse powers of the high scale A: 

J? = J? SM + 5^ + 5^ e + --- , with JjBioc^O'. (3) 

After EWSB, some of the effective operators may result in corrections to the low-energy 
SM parameters. In addition, new exotic couplings may result, such as those in Eq. (1). 

The only possible d = 5 operator is the famous Weinberg operator [6], which leads, after 
EWSB, to Majorana neutrino masses. We will not need to consider it for the present study 
of NSI. 

At d — 6, some operators modify the low-energy standard couplings, among which non- 
unitary corrections to the leptonic PMNS mixing matrix are specially relevant to the case 
under study. Departures from unitarity are a general feature of beyond the SM models 
involving exotic fermions [2,7-15]. All fermions of the same charge will generically mix 
through the mass matrix, leading in those theories to unitary mixing matrices of dimension 
larger than three, while the effective 3x3 sub-matrix - relevant at low energies - is not 
unitary. In the effective Lagrangian formalism, the effect appears technically at leading 
order through d = 6 gauge invariant operators involving only two fermions, which induce 
non-canonical corrections to the fermion kinetic terms [7-11]. Such operators are for in- 
stance typical of fermion-mediated Seesaw scenarios. The trademark of non-unitarity is 
that the coefficients of the NSI operators induced by it and contributing to neutrino pro- 
duction, detection, and matter effects are not independent but related. Barring fine-tuned 
cancellations, the stringent bounds and future signals on non-unitarity [16-18] apply as 
well to NSI, except for those NSI operators affecting exclusively the propagation in matter. 
Recently, the value of the elements of the PMNS matrix have been extracted from data 
without assuming a unitary mixing matrix [19], and new related CP-odd signals have been 
proposed as well [20-22]. For a detailed discussion of the NSI-non-unitarity relationship, 
see [23] . We will leave non-unitarity out of the main line of this work, discussing only its 
qualitative implications. 

Effective interactions such as the NSI in Eq. (1), obviously require to consider operators 
made out of four leptonic fields, plus Higgs fields in the case of operators with d > 6 [24-26]. 
There is a plethora of d = 6 [27] and d = 8 [28] operators 1 , with different classes of models 
resulting in different operators and operator coefficients. Those among them relevant for 
NSI can affect neutrino production or detection processes, or modify the matter effects 
in their propagation, depending on the operator or combination of operators considered. 
Notice, though, that the coefficient of a d — 8 operator is expected to be suppressed by a 

d = 7,9... operators [29] are odd under baryon and lepton number and not relevant for the present 
discussion. 
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factor v 2 /A 2 with respect to d = 6 operator coefficients (where v is the vacuum expectation 
value of the Higgs field (H°) = v/y2~ = 174 GeV), and thus negligible unless the new scale 
is very close to the electroweak one. 

Much effort has been dedicated to analyze the experimental constraints and future detection 
prospects of NSI [24,30,31]. This encompasses their impact on weak decays, solar and atmo- 
spheric neutrino physics [32-39], astrophysics [40,41], early universe processes [42], collider 
and neutrino scattering experiments [28,43-46], and past and future neutrino oscillation 
experiments including neutrino factories [47-69]. Overall, the signals involving the v T field 
are at present the least constrained ones. 

We emphasize that we study "large" NSI, where "large" refers to being potentially observ- 
able at future neutrino oscillation experiments. For example, for the flavor- changing NSI 
interfering with standard oscillations, bounds 10~ 2 to 10™ 3 (relative to the SM four-fermion 
interaction coupling) on the operator in Eq. (1) can be expected from a neutrino factory 
(see, e.g., Ref. [65]). After EWSB, the four-fermion interactions from Eq. (3) will be sup- 
pressed by roughly {v/ K) d ~ i with respect to the SM four-fermion interactions, which means 
that A is allowed to be largest for the d = 6 effective operators, whereas for d — 8 oper- 
ators, it has to be very close to the electroweak scale in order to produce a sizable effect. 
Therefore, we focus on the d = 6 operators first, and then increase the level of complexity. 
However, from this simple comparison, one can already read off that new physics above the 
TeV scale will be very difficult to be observed at future neutrino oscillation experiments, 
since the suppression with respect to the SM is roughly 10~ 2 (for d — 6) and 10 -4 (for 
d — 8), respectively - especially if the d = 6 operators turn out to be not good enough. 

Most of the literature deals mainly with NSI in experiments, and does not discuss the mod- 
els behind them - with some exceptions [3,5,23,25,26,70-73]. In this work, we will classify 
all d = 6 and d — 8, SU(3) x SU(2) x U(l) invariant, leptonic NSI operators in terms 
of the heavy mediators inducing them. This can be achieved in a model independent way 
decomposing each operator into all possible products of currents. The SM charges of the 
corresponding field combinations will be the SM charges of the putative heavy mediators. 
The study will thus be confined to the case in which the tower of non-renormalizable opera- 
tors has been produced by the tree-level exchange of new heavy fields, whose precise nature 
other than their SM charges we do not need to know. 

We will first emphasize the case in which the mediators exchanged only couple to SM bilinear 
field combinations. In this study, we refer to "SM bilinears" as fundamental interactions of 
exactly two SM fields with one or two exotic fields, where the latter possibility amounts to 
couplings between two exotic bosons and two Higgs doublets. Other than that, there can be 
in addition new exotic couplings involving only one SM field, which will also be addressed 
later on. 

The decomposition in terms of SM bilinears, initiated in [25,26], immediately leads to 
correlations between previously uncorrelated effective NSI operator coefficients, with very 
fruitful physics consequences, as we will show below. For instance, it has been realized 
that the lowest dimension operators, that lead to NSI without simultaneously inducing 
dangerous transitions among four charged leptons, are d = 8 operators. The first example 
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proposed [24, 74] is of the form 

NSL = (L i H i )Y{H^ i L i ){Ej p E), 



(4) 



where i, L E and H denote the SU{2) L index, leptonic doublet, leptonic singlet and Higgs 
doublet, respectively (note that we use the convention with Y — — 1/2 for the Higgs doublet, 
H T = (^H°, H~)). We will illustrate below, though, that any realistic model responsible 
for it will also induce other dangerous d = 6 operators and/or some other low-energy 
effects for which stringent experimental bounds may exist (i.e. non-unitary corrections to 
the PMNS matrix, corrections to the electroweak precision data, flavour- changing neutral 
currents (FCNC), etc.). Eq. (4) is only one example of the NSI seeked. We will determine in 
this work several other independent d — 8 operators which lead to large NSI and no charged 
lepton processes. The symbol Onsi will be extended to denote generically any element in 
this ensemble. 

In our general analysis, after determining all possible mediators, the resulting correlations 
between the possible d = 6 and d — 8 operators will be systematically studied. We will then 
establish which mediators or combinations of mediators can lead to large NSI, without in- 
ducing experimentally excluded leptonic charged flavour- changing transitions, and/or other 
undesired phenomenological consequences. 

Our main motivation in this study is to determine the minimum level of complexity needed 
for a viable model of NSI. As an illustration for model building, a particular simple toy 
model will be developed in which the operator Cnsi above is induced unaccompanied by 
any leptonic d = 6 operator. The aim is to show the generic prize to pay at the theoretical 
level for allowing observable NSI effects at future experiments. 

Note that we focus in this study on the necessary conditions to build a model with large 
NSI, while for any given model additional limitations may arise. Supplementary constraints 
which could arise from a phenomenological analysis at one-loop are also not considered here 
and should be addressed when considering a particular model. From the experimental point 
of view, we will not make any explicit statement how likely it is to observe large NSI. We 
leave the interpretation of this likeliness by judging the necessary conditions for a viable 
model to the reader. Finally, possible NSI involving quark fields are neither included in this 
study. 

2 Effective operator formalism 

The SM Lagrangian is extended to accommodate the tower of effective operators 



where the two terms run over all possible d = 6 and d = 8 operators relevant for purely 
leptonic NSI. The flavour composition will be made explicit in each coefficient and operator, 
i.e., {Ci)°p S (Oi)^L. All distinct flavour combinations for the same operator structure will be 
taken into account, as they correspond in fact to independent operators. 




(5) 
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Effective operator basis 



In order to find all possible d — 6 and d = 8 effective operators leading to purely leptonic 
NSI, we will use the following bases: 

• d = 6 operators. A complete basis of d = 6 operators invariant under the SM gauge 
group and made out of the SM light fields was proposed by Buchmuller and Wyler 
(BW) [27]. The four fermion operator structures relevant to our problem are: 



(O^H^E^E 3 ^), (6) 

(Ol^HL'YL^L^L,), (7) 

=(L P YrL a )(L s lp rL,) , (8) 

{0 EE )^={E^E a ){E s lp E,), (9) 



where L (E) refers to the SU(2) leptonic doublet (singlet). We will refer to the 
coefficient matrices for these operators by (C LE )p], (Cl L )p], (C| L )^, and (C EE )p], 
respectively. The operators O ee do not produce NSI directly, but will play a role 
when considering charged lepton flavor violation (since they share some mediators 
with the operators in Eq. (6)-(8)). 

On top of the above, there are two d = 6 operator structures including two lepton 
doublets L and two Higgs doublets H, 

{0\ H t = {L^H)\p{H ] L a \ (10) 
{Ol H )i = {L^H)iip{H^L a \ (11) 

and a operator with two £"s and two if's 

{0 EH )i = (HhD<>H) {E%E a ) , (12) 

where D denotes the SM covariant derivative. These three operators belong to the 
class which, after EWSB, correct the parameters of the SM Lagrangian. In particular, 
they renormalize the kinetic energy of neutrinos and/or charged leptons [9-11]. As 
previously mentioned, they result in non-unitary corrections to the the leptonic mixing 
matrix and/or correct the charged and neutral electroweak currents [27], and will 
not be further developed in this work. We include them above only for the sake of 
completeness (see also Sec. 4. 2. 2). 

• d = 8 operators. A basis was discussed by Berezhiani and Rossi (BR) [28], with the 



relevant operators given by 

{Ol EH )% HL P l P L a ){E s lp E,) (rfH) , (13) 

{0*leh)% =(L"YrL a )(E\E,) {H^tH) , (14) 

(Oll^ =(L P YL a )(L s lp L,)(H^H), (15) 



5 



(02&r)£ =(L"Y?L a )(L\rL y )(H^H), (16) 

(C>iEr)£ ={L P YL a ){L & lp fL,){H^H\ (17) 

(OE5r)£ ^iV^KiV,)^), (18) 

(Of^)^=(-ie o6c )(Z^ 7 W a )(ZV 6 ^)(^ t ^) ) (19) 

(O eeh )^ =(EYE)(E lp E)(H^H) . (20) 



In these operators, subscripts correspond to a shortcut notation for their SM field 
composition, whereas superscripts denote the corresponding £77(2) charges of the 
field combinations. Once again, although the operators Oeeh cannot induce NSI by 
themselves, they will come to play a related role, as they induce charged lepton flavour 
violating transitions. 

Strictly speaking, not all of the above operators are independent when the full flavor 
structure is taken into account, as 

{Ol? H fl = {011^1. (21) 

However, the expressions below will look much simpler if both operators are used. 

Notice that the phenomenologically interesting (9nsi operator in Eq. (4) can be ex- 
pressed as a combination of the two first operators in the list above, 

0nsi = \ (Ol EH + Ol EH ) . (22) 

This means for instance that if a model only induces at d = 8 the operators 0\ EH 
and O leh with similar weights and no d = 6 NSI operator, it could be an optimal 
candidate for viable large NSI. We will explore later some examples of this kind. 



Decomposition in terms of SU(2) field components 

After EWSB, the contributions from the d = 6 and d — 8 gauge invariant operators result in 
two very simple sets of operators. From the LLEE-type operators, Eqs. (6), (13) and (14), 
we find: 

1/1 v 2 \° 7 

6J? eS =— (--C LE + ^(Cl EH + CIeh))^ (*VPz^) ?7pP^ 7 ) 

+ Y 2 [~2 Cle + ^ {Q1leh ~ Cleh) ) 8 (^ P ^) (^pPiA) + h - c - » ( 23 ) 



The first line in this equation produces the relevant NSI, whereas the second line leads to 
the (unwanted) four charged lepton contributions. The NSI in the first line involve only 
right-handed charged leptons. In consequence, their effect at the neutrino source will be 
chirally suppressed 2 . 



2 At detection, the effect of these NSI is subdominant because of the dominance of the neutrino- nuclcon 
cross section. 
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From the operators involving four lepton doublets, Eqs. (7), (8), (15)-(19), it results that 3 
^ cff = - (« LL )^ (*VP^«) (^7 p Pl^ 7 ) 

+ -1 (cl + ci L + ^ (ci- + cf- - ca - cs&))^ prPx^) (ivpza) 

+ h.c. , (24) 



where 



/ 2 \ a / 

nLLLL\ al _ I n\ _ r>3 , _^ ( n\W _ r 331 _ A" 133 i (^313 \ | 

U NSI J p S ~ y^LL u LL~r~ 2y^2 \ LLH U LLH U LLH ~r ^LLHJ J 

, { n 1 r 3 i ^111 ^-331 , ^133 ^,313 

"r 1 U LL U LL "r ^2 \ LLH U LLH "> U LLH U LLH ) I 

i ^r 3 + ^!_(? 7331 _r 333 ^ 



S/3 



Sf3 

A^ 1 



+ [2Cl L+ V -{Clf H + ClT H )\ . (25) 



,36 



Note the different flavor structure in the four lines in Eq. (25). In addition, note that the 
term relevant for the NSI, i.e., the first line in Eq. (24), couples to left-handed charged 
leptons, which means that source NSI can be generated as well. In resume, matter NSI are 
(not) correlated with source and production NSI for LLLL (LLEE)-type operators. 



Connection to NSI and phenomenology 

Let us first consider NSI in matter. The phenomenology of neutrino propagation under these 
conditions is customarily described in terms of the Hamiltonian in the flavour basis [4,24,51], 

if/ \ A + e ee C M e£\ } 

H F = ^{U\ Ami & + «cc (e£ )* e™ & ]\, (26) 

I V Amy V it™)* KT J 

where, acc is the usual matter effect term defined as acc = 2\Z2EGfN £ (with N e the 
electron number density in Earth matter). 

From Eqs. (23) and (24) it follows that 

2 

m,L _ v ( n LLLL\ ae m,R _ v L n v / n l , n S 



~p a 2A 2 



{ C ^)p e ' ^ = 2V \2 Cle + 2A^ (C ^ + CleH % £ ' (2T) 



3 Here we do not show the interactions among four neutrinos which these operators also induce. See 
Appendix A for a discussion of these interactions. 
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with C^[g I ii as defined in Eq. (25). These two parameters in matter lead to a total 

jm rn.L . m,R /OQ\ 

because matter effects are only sensitive to the vector component. 

In addition to the propagation in matter, the production or detection processes can be 
affected by NSI. For the specific case of a neutrino factory and considering just the purely 
leptonic NSI under discussion, only effects at the source are relevant, since the detection 
interactions involve quarks 4 . They are customarily parametrized in terms of e s a/3 , which 
describes an effective source state \u^) as [4,48,75] 

K) = K)+ E e ^K>- ( 29 ) 

In this case, the muon decay rate could be modified by the NSI interaction in Eq. (24), with 
the largest effect resulting from the coherent contribution to the state at the source [48,53]. 
It appears as an admixture of a given flavour v a with all other flavours, encoded by v 1 in 
Eq. (29). Two types of contributions are possible, 

2 2 
c s _ (pLLLL\e^ ™~ c s _ _^ ( nLLLL\^e /on\ 

M/3 ~~ 2A 2 P e 13 ~ 2A 2 1 13(1 ' 

The second possibility will affect the golden v e — > appearance channel, where the effect 
might be easiest to observe. If the coefficients in Eq. (25) are known for a specific model, 
one can easily calculate the connection between source and propagation effects via Eqs. (27) 
and (30), a connection which does not hold for LLEE-type operators above, as explained 
earlier. 



Conditions to suppress charged lepton processes 

Let us discuss now potentially dangerous contributions to charged lepton flavour violation 
processes, possible modifications of Gp and the constraints on lepton universality. The 
focus is set on pure charged lepton processes at tree level. These interactions can result 
from the second terms in Eqs. (23) and (24). They should be very suppressed in any 
phenomenologically viable model. In order to cancel those terms, the putative beyond the 
SM theory has to satisfy, to a high degree of accuracy, the following constraints: 

/ 1 v 2 \° 7 

\2 Cle + Yk^ cXhEH ~ c1eh) ) 5 = ' (31) 

(2 \ Q 7 

Ch +CL + ^- 2 (Cl? H + Clf H - C\f H - CH%) ) r 0, (32) 

for all possible values of the flavour indices (Greek letters). A possibility suggested by these 
equations is that there could be cancellations among d = 6 and d = 8 operator coefficients. 

4 Superbeams, for instance, use hadronic interactions for neutrino production, which are not affected by 
purely leptonic NSI to first order. 
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However, we will not discuss such a possibility in this study, as it would correspond to fine- 
tune the scale A. We will therefore require that the d = 6 and d = 8 operator coefficients 
in Eq. (31) and Eq. (32) cancel independently. 

For the d — 6 operator coefficients, it reads (omitting flavor indices) 

Cle = , C LL = —Cf L , , C E e = , (33) 

which implies that only LLLL-type operators can induce large NSI. One possibility for 
its implementation is the antisymmetric operator mediated by a SU(2) singlet scalar in 
Ref. [26], which turns out to be the only d = 6 possibility requiring just one tree- level 
mediator, as we shall explicitely demonstrate. 

For the d — 8 operator coefficients, the cancellation conditions read 

C-leh — Cleh > Cllh + C-Tlh - ClTh - Cllh — > C-Tlh arbitr. , C E eh = , (34) 

where the first condition corresponds to operators of the type Cnsi m Eq. (4), see Eq. (22). 
In the following, we will refer to operators satisfying Eq. (34) as Onsi ; &.e., we define the 
class of potential non-standard neutrino interaction operators in mass dimension eight as 
the one which does not introduce any harmful d = 8 processes with four charged leptons. 
Eq. (4) is (apart from Fierz rearrangements) the only such possibility with two right-handed 
charged leptons involved. When considering leptonic NSI involving four left-handed fields, 
several new operators of this kind will be determined later on. 

As far as the possible NSI in terms of SU (2) field components are concerned, not all flavour 
structures can be generated from the d = 6 effective gauge invariant operators if charged 
lepton processes are suppressed. Applying the d — 6 cancellation conditions in Eq. (33) to 
Eq. (25), it results that the d = 6 contribution to the coefficient (^s/" L )^J is antisymmetric 
in the flavor index exchanges (0,7) — > (7,0) and ((3,5) — > (5,(3), which means that a 7^ 7 
and (3 7^ 5 for viable NSI. As regards matter effects, this implies that only e™^, e™, and 
- defined in Eq. (27) - can be generated from d = 6 operators and the connection with the 
source effects is given by 

e A*M = ~ e ee = ~ e (J,fJ, ) (35) 
= -«r)* • (36) 

In contrast, e™ is not connected to the source effects at the effective operator level 5 . Notice 
that, for instance, the NSI in Eq. (35) contribute to the Gp measurement coherently (i.e., the 
interference with SM couplings contributes linearly to the rates), for which quite stringent 
bounds exist. These results hold in general for any purely leptonic NSI d = 6 operator with 
suppressed interactions among four charged leptons, i.e., Eq. (33). Furthermore, for the 
particular case of a neutrino factory, the antisymmetry conditions and constraints described 
above imply that the only possible non-negligible NSI source terms induced by d = 6 
operators are e s eT and e^ T . 



5 There can be also subdominant effects in detection chains. For example, in OPERA, the taus resulting 
from hadronic interactions decay into muons or electrons. It implies for instance e^. = — e s TT1 which means 
that tau decay into electrons is in this case connected with matter NSI. Note that NSI and SM contributions 
add coherently to the r~ — > e~v T v e width. 
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A T GTT{<1\ 

l\L bU{Z)L 


U(1)y 


Lorentz 


Mediator 


Bilinear(s) 


Models [fteis.J 


2 1 


— 1 


scalar 


Hi 


L c \t L 


Zee model [76-78], ^pSUbY [79] 




-2 


scalar 


I* 2 


We 




3 


-1 


scalar 


3ii 


T?\T 2 T a L 


Left-right sym. [80-83] 


2 


-3/2 


vector 


Z -3/2 




331 model [84-86] 


1 





vector 


-IV 

x 


LjPL, EjPE 


Models with Z' [31] 


3 





vector 


Q V 

d 


LYr a L 


Models with W' [31] 


2 


1/2 


scalar 


Z l/2 


EL 


%SUSY [79] 



Table 1: Possible SM bilinear field combinations involving only leptons. AL refers to the lepton number 
of the mediator, SU(2)l to electroweak singlets (1), doublets (2), or triplets (3), U(l)y to the hypercharge 
Y = Q — 1^ , and "Lorentz" to the Lorentz nature. The mediator notation is Xy, where X, C and Y 
denote its SU(2), Lorentz, and Y properties, respectively. The mediators which carry two units of lepton 
number were studied in Ref. [87]. 

3 Model analysis of d = 6 operators 

In this section, we discuss the model-building implications of requesting large d — 6 NSI 
induced by theories of physics beyond the Standard Model. We specifically highlight the 
basic principles, which can be found in the d = 8 well. However, as we shall see 

later, the d = 8 case is technically somewhat more challenging. 

In order to shed light on model building, let us analyze the operators according to the 
possible tree-level mediator. This is most efficiently done by listing all possible SM bilinear 
field combinations, and combining them in all possible ways [25,26]. 

We therefore show in Table 1 the possible bilinears constructed from leptons only, which 
can lead to the d = 6 NSI operators in Eqs. (6)- (9). It is obvious from the table that the 
bilinears carry the mediator information and that they can therefore be directly associated 
with specific models (as illustrated). The mediators are denoted - all through the paper - 
by Xy, where 

• X denotes the 577(2) nature, i.e., singlet 1, doublet 2, or triplet 3. 

• C refers to the Lorentz nature, i.e., scalar (s), vector (v), left-handed (L) or right- 
handed (R) fermion 6 . 

• Y refers to the hypercharge Y = Q — 1^ . 

Table 2 shows in turn all possible d = 6 operators which can be constructed from the SM 
bilinear field combinations in Table 1. The coefficients of the d = 6 operators obtained 
by this procedure are denoted by (c X£ )^J and (/ X£ )^J, where c (/) indicates that the 
corresponding operator results from the exchange of particles carrying two (zero) lepton 

6 Fermionic mediators will appear explicitly later on, when discussion d — 8 effective interactions 
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d — 6 operators Mediator Cle C\ l C\ l Cee 

LEEL 



(c 2v /A 2 )((E c )^p L„) ( ZAy„( E c Y) 

(flv 1 \2\(f P-pT UffL rp \ 

(f 2s /A 2 )(^E,)(E s L a ) 


-3/2 
1 11 

J-o 

Z l/2 


2c 2v 
~ z Jle 






L L Tj L 










(c 3 VA 2 )((L c ) a zr 2 fL 7 )(L' 3 fzr 2 (L c ) 5 ) 

(fll/A 2 )(L^YL a )(L\L 7 ) 

(f^/A 2 )(L^rL a )(L s lp rL,) 


Hi 

3ii 

°0 




4 L 
3 „3s 
4° 
flv 

Jll 


l r ls 
4 C 
l„3s 
4 C 


EEEE 










(c^/A 2 )((^)^ 7 )(^(^ c ) 5 ) 
(/M(BV4)(£Vt) 


1 s 






fit; 
J EE 



Table 2: Possible d = 6 operators obtained by combining the lepton bilinears in Table 1. The coefficients 
are labeled c X£ (AL = 2) or / X£ (AL = 0), where the X and C are defined as in Table 1 and the subscripts 
refer to the combination of bilinears, in an obvious notation. The last four columns contain their contribution 
to the d — 6 operator coefficients in the BW basis in Eqs. (6)-(9). The flavour structure for any coefficient 
in the table is understood to be ()jgj, see main text for further explanations. 

number, and X and £ refer again to the SU(2) and Lorentz nature, respectively. Any 
subscript refers to the combination of bilinears involved. 

At this point it is important to note that the operators obtained from the mediators do not 
constitute a basis. Instead they are not independent, but linear combinations of those in 
the BW basis, Eqs. (6)-(9). Therefore, it might be more accurate to call them "mediator- 
operators" or "operator combinations". We will not make this special distinction, but the 
reader should keep that in mind. Re-writing the individual effective operators from Table 2 
in the BW basis, we find the coefficients given in the last four columns of Table 2. For 
example, the first line of the second group, mediated by l s _i, reads (including flavor indices) 

( r ls \»y 1 ( Is scry ( Is \<*7 

^ef 6 = m^L^T^f) = 1^4* (QUI* _ 1 MfM. {0 l L )^ . 

^ (37) 

Conversely, the decomposition of the operator 0\ L of the BW basis can be read off from 
the column labeled C\ L , in terms of the relative weights of the mediator-operators: 

{Cl L )f 5 = \{<&)% - \{c*%] + (fH)«] . (38) 

Note that the flavor indices in the first column of Table 2 are arranged such that the flavor 
indices of all coefficients and of the BW operators are the same as in Eqs. (37) and (38). 
Therefore, we show the flavor indices explicitly only in the first column. 

In order to have large NSI without four charged lepton interactions, the d = 6 cancellation 
conditions Eq. (33) must now be implemented. One can directly read off now from Table 2 
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1/Ml j 1/Ml ] 1/Ml 

Figure 1: Diagrams mediated by a bilepton 11 x . The effective dimension six operator results from the 
the first diagram at energies below the mediator mass M± s . However, the fundamental interaction will also 
create the diagrams in the middle and right, and the couplings will be related. 

that these conditions can be re-written as 

2c 2v - 2fH + f 2s = (from C LE = 0) , (39) 
~c 3s + fll + f 3v = (from C\ L + C\ L = 0) , (40) 
Cee + Vee = (from C EE = 0) , (41) 

in the mediator picture. The operators contributing to the first equation will not produce 
any NSI (since Cle = in Eq. (23)), while the operators present in the second equation 
lead to NSI if C\ L = -C\ L £ (c/., Eq. (24)). 

One approach to use Table 2 is to discuss departures from the SM couplings. For example, 
for a hypothetical experimental departure pointing towards a four-lepton coupling such as 
that in operator 0\ L in Eq. (8), Table 2 indicates directly that a new heavy scalar triplet 
could induce it at tree-level, while a scalar doublet wouldn't. 

From the model building perspective, it is illustrative to consider again the case of the 
operator mediated by l s _ 1 leading to Eq. (37). The table shows that it is the only d = 6 
possibility using only one mediator which directly satisfies the cancellation condition of pure 
charged lepton interactions Eq. (33) (or their tree-level equivalent Eqs. (39)-(41) ). This 
antisymmetric combination of the basis elements was first found in Ref. [26]. 

This example serves to illuminate the power of the mediator analysis (see Ref. [26] and 
also Ref. [23]). The l s _ 1 exchange leading to the originally proposed operator is depicted 
in Fig. 1, left. Once a certain mediator is assumed for a certain operator, contributions to 
other operators are simultaneously induced, though, as illustrated at the center and right 
of Fig. 1, i.e., 



( C il)erl _ I^W| 2 | A er | 2 _ | ( c il)e^| l( C il)erl 



A 4 M{ s A 4 



(42) 



where \ a p is the coupling for the lepton-bilepton interaction and M± s is the mass of the 
bilepton. A coherent contribution to Gp and a violation of the lepton universality is then 
induced by the diagrams at the middle and right of Fig. 1. From the strict experimental 
bounds on these quantities has been constrained to |e™ | < 1.9 ■ 1CT 3 (90% CL), using 
this particular mediator [23]. The bound from a neutrino factory on |ej^_| would be 1.8 • 10~ 2 
for complex 6™ [65]. If it was assumed to be real, which does not describe the most general 



12 



class of models, the bound would be 3.7 • 10 4 (90% CL) [65]. However, since this is a 
model-dependent assumption, we do not use this bound. 

The antisymmetric operator discussed in the previous paragraphs is not, however, the only 
possibility to build a model satisfying Eqs. (39)-(41). For example, one may choose bosonic 
triplets 31 1 and 3q, for which the coefficients can be chosen independently, in order to satisfy 
Eq. (40) without suppressing completely the d = 6 NSI operator coefficient. In particular, if 
the simplest possibility is experimentally constrained, one may consider models with more 
than one mediator. 

At this point, we would like to clarify that cancellations or fine-tuning of operator coefficients 
cannot be an argument in itself for judging the naturalness and complexity of a model. It 
depends on the field content and the symmetries of the model. Consider for instance once 
again the antisymmetric operator in the left-hand side of Eq. (37), induced at tree-level by 
the exchange of just one mediator, 11 1; illustrated in Fig. 1 left. That equation shows that, in 
the BW basis, the antisymmetric operator appears to be constructed from the combination 
of two BW operators with specific (fine-tuned?) coefficients. In the effective operator 
picture, "fine-tuning" is thus a basis and model-dependent qualification. We therefore 
define the simplest model to be the one with the fewest mediators. In the d = 6 case, it is 
the antisymmetric operator in Eq. (37) with only one mediator. In the case that the NSI 
come only through d — 8 (or higher dimension) effective operators, we will demonstrate 
that the simplest viable models require at least two mediators. Once the field content is 
chosen, any relative precise adjustment of the couplings of the mediators can be considered 
a fine-tuning, unless the symmetries of the model ensure it. It will be left to the model 
builder to eventually explore possible symmetries, whenever such cancellations will turn out 
to be required below for phenomenologically viable NSI. 

In general, it is easy to show that all NSI from d = 6 operators are strongly constrained when 
the possible mediators are taken into account. There is, however, one exception. The present 
experimental constraints allow the condition Cle = in Eq. (33), which cancels interactions 
among four charged leptons, to be substantially violated for certain combinations of flavour 
indices. In particular, the coefficient of the flavor conserving (BW) operator (L T E e )(E e L T ), 
is not very strongly constrained [28,74]. The mediators 2 v _ 3 / 2 or 2*/ 2 (c/., Table 2) can 
generate such an operator, leading to the following effective interactions, cf., Eq. (23): 

ssg g* = ((pT-fP L v T ) (e lp P R e) - (fYP L r) (e lp P R e)) + h.c. . (43) 

The coefficient is constrained by (see Eq.(14) in Ref. [28]) 

\e™\ = ^ 2 \{C LE y T l\ = \K TR \<Q.l. (44) 

If the possibility of large SU(2)i breaking effects was considered in addition, a possible 
gain of almost an order of magnitude could be allowed for the NSI e™ strength [26]. In 
conclusion, large (order unity) values for e™ resulting from d — 6 effective interactions are 
not excluded. 

Table 2 also shows that the relationship between mediator and coefficient is unique at the 
d = 6 level, except for 1q. If a model uses this mediator, then there will be three different 
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d = 6 operator contributions, which are independent in the BW basis. In particular, one 
cannot neglect Oee, which can induce physics effects while not resulting in NSI. 

4 Model analysis of d = 8 operators 

We consider all possible d = 8 operators which can induce purely leptonic NSI, analyzing 
them from the point of view of their possible tree-level mediators. We will focus on the sys- 
tematic analysis of all possible products of SM bilinears, which may result from exchanging 
mediators which only couple to pairs of SM fields. Such an analysis was performed for d = 6 
operators in Refs. [25,26], and we extend it here to the d = 8 case. Other scenarios leading 
to some of the d = 8 operators will be briefly analyzed afterwards. 

A convenient basis of linearly independent d = 8 operators has been given in Eqs. (13) 
to (20), i.e., the BR basis. In order to suppress four charged lepton interactions, both the 
cancellation conditions for d = 8 operators in Eq. (34) and the cancellation conditions for 
d — 6 operators in Eq. (33) are now required to be satisfied. Under these conditions, if 
any d = 6 operator remains, it is expected to dominate the new physics and, as discussed 
in the previous section, only effects related to e™ are then allowed to be experimentally 
sizeable. In this section, we instead focus on NSI which stem exclusively from d = 8 (and 
higher) operators and their implications for model building. In particular, we are interested 
in the Onsi without four charged lepton interactions, i.e., satisfying Eq. (34), which has 
been object of intense speculations in the literature. 

When the mediators couple only to SM bilinears we have the following options with respect 
to the undesired d = 6 operators: 

1. The required mediators do not induce any d = 6 operator involving four leptons (in 
other words, the mediators differ from those in Table 2). 

2. The d — 6 couplings induced by different mediators turn out to explicitly cancel among 
themselves. 

As we will illustrate later, there is no simple possibility for which the first option works. 
For the second option to happen, the coefficients for the BW operators in Eqs. (6)-(9) have 
to vanish independently, because they constitute a basis: 

C LE = 0, Cl L = 0, C 3 LL = 0, C EE = 0. (45) 

Their implementation in the mediator picture can be read off from the columns in Table 2. 
They are given by Eqs. (39) and (41), together with 

\c\ s l ~ ~/ s + fll = (from C\ L = 0) , (46) 

-\<&-\*' + f*' = (from C» L = 0) , (47) 

which replace Eq. (40) of that set. For example, if a model introduces two bosonic doublets 
2^ 3 / 2 and 2* , 2 , one can satisfy Eq. (39) (to which Eq. (45) simplifies in this case) by achieving 
2c 2 " + f 2s = 0. 
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Note that the introduction of exotic fermions in the game potentially leads to the additional 
d = 6 operators in Eqs. (10)-(12), which are made out of two lepton fields and two Higgs 
doublets. In accordance with the main line of this section, we do not consider constraints 
from those operators, which means that, unless explicitly stated otherwise, when mentioning 
d — 6 operators in this section we refer exclusively to those in Eqs. (6)-(9). 

4.1 A toy model 

In order to estimate the theoretical price to pay for obtaining large NSI from exotic particles 
coupling to SM bilinears, without large charged lepton flavour violation, we show here a toy 
model in a bottom-up fashion, which precisely generates the d = 8 operator Onsi i n Eq. (4) 
and no d = 6 operator. Then we will provide a systematic analysis, from which we will 
recover the toy model as the simplest possibility in a top-down approach. 

Consider the following toy Lagrangian for the underlying theory, which adds both a new 
scalar doublet (2^ 2 ) $ and a vector doublet {2 V _ 3 ^ 2 ) to the SM Lagrangian, with general 
couplings to the SM fields y, g and A's, 

2 = ^sm - (y)^ (1?)%^ - (g) ps {L p ) i Y{E c ) s {V p )i 

+ \ lv {H ] H){VlV p ) + \ 3v (hWh)(V}tV p ) + h.c. + ... (48) 

where the dots refer to other bosonic interactions not relevant for this work. After inte- 
grating out the intermediate particles, the following d = 6 effective interactions involving 
leptons are induced (see Table 2): 

^eff =6 = (E^YL^L^E^) + (L?E 7 )(E s L a ) , (49) 

where now 

A 2 My 'A 2 M| { } 

For simplicity of notation and illustrative purposes we can assume Mq, ~ My = M(= A). 
The d = 6 cancellation conditions on four charged lepton transitions in Eq. (33), or its 
equivalent in the mediator picture Eq. (39), translate into 

-2(gty<*(g)p S + (yt) s a ( y y = 0. (51) 
The relevant effective d = 8 Lagrangian induced reads 

^cff =8 ~[Ai, (LyE) (Ey^L) (&H) + \ 3s (LyE)r(Ey^L) (H^tH) 

+X lv (LgYE c ) (E c l P g ] L) (H*H) + A 3 „ (LgYE c ) r (W lp g^L) (H*tH)], (52) 

where flavour indices have been omitted and each expression in brackets is to be understood 
as a flavour singlet. Eq. (52) can be rewritten in terms of the operators of the BR basis in 
Eqs. (13) and (14), as 

d-^eS 8 = _ ]\4 (^LEH^LEH + ^LEH^LEh) j (53) 
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Figure 2: Dimension eight operator decomposed into dimension four interactions 



where 



C 1 



c 3 



LEH 



A 3 , (g ] V a (g)ps + \\3 S (y j )s a (y)i' 



(54) 
(55) 



In order not to produce interactions between four charged leptons, it is necessary to satisfy 
Eq. (34), i.e., the condition C\ EH = C LEH ^ 0, so that the effective d — 8 interaction 
in Eq. (53) reduces precisely to (9nsi in Eq. (22). Together with the d = 6 cancellation 
condition, Eq. (51), it is finally required that 



(56) 



As a consequence, the NSI in matter can be substantial for all flavours. While source and 
detection NSI cannot be created from our toy model, the epsilon matter parameter reads 



m,R 



2M 4 



\(X ls + X lv )(g^r(g)p, 



(57) 



In resume, by adding both an SU (2) doublet scalar and a doublet vector to the SM content, 
and imposing two relations to their couplings, Eqs. (51) and (56), a toy model for viable 
large NSI has resulted. The model interactions are visualized in Fig. 2, where the first two 
effective interactions in Eq. (52) correspond to the diagram on the left - mediated by 2\ , 2 
- and the last two interactions to the diagram on the right - mediated by 2 v _ 3 / 2 - m fact, 
other combinations of just one of the first two operators in Eq. (52) together with one of 
the last two operators in that equation would have been enough for the purpose 1 . As we 
will demonstrate below, our toy model is the most general possible model involving only 
two mediators, when the exotic particles couple only to SM bilinears. 

We keep dubbing the construction above as "toy" because, to begin with, the presence of 
a vector field, which is not a gauge boson, implies that it is non-renormalizable. The toy 



7 For instance a combination involving \\ s and A3t,, or alternatively A3 S and Ai^, would be suitable. 
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Lagrangian, Eq. (48), can thus only be considered as an effective theory of some larger 
construction, such as for instance models of extra dimensions in which the vector doublet 
could be a component of a higher dimensional gauge theory. 

Moreover, its phenomenological analysis is beyond the scope of the present work: the con- 
straints from electroweak precision tests need to be analyzed for each specific model, in 
particular the oblique corrections [88-90] it may induce. The new couplings may also have a 
relevant impact on other flavour changing transitions at the loop level, although considering 
large values for the quartic couplings A and small values for the elements of the g and y 
flavour matrices, it will probably remain phenomenologically safe. 

The toy model demonstrates that it is possible to achieve the desired d = 8 interactions, 
without inducing simultaneously d — 6 ones, by fixing the coefficients of the new fields in 
the Lagrangian. It requires ad-hoc cancellations, though, and it is left as an open question 
for the model builder whether some symmetry can justify them. 

4.2 Systematic analysis 

In this subsection, a systematic analysis of all possible effective NSI d = 8 operators is 
performed. The full decomposition of any combination of d = 8 operators, constructed 
from combining bilinear combinations of SM fields, leads to a large number of possibilities. 
We will first consider the cases which are conceptually similar to the toy model above, i.e., 
new fundamental interactions involving exactly two SM fields, which are the SM bilinears 
according to our earlier definition. Then we will discuss new interactions involving only one 
SM field. 

4.2.1 Mediators coupling to SM bilinears 

We summarize these possibilities for the LLEE-type operators in Table 3 and for the LLLL- 
type operators in Table 4, which are one of the main results of this study. The notation 
used has been described in Sect. 3. The tables show, from left to right in each row: 

• an ordinal assigned to each operator, 

• the operator itself, 

• the value of the operator coefficients of the BR basis needed to reconstruct it, 

• whether the d = 8 cancellation conditions in Eq. (34) are directly fulfilled ("Cnsi?"); 

• the required mediators, with those inducing additional d = 6 interactions of four 
charged leptons (Table 2) highlighted in boldface. 

Obviously, the number of possible mediators of d — 8 interactions is much larger than for 
the d = 6 case in Table 2. In particular, fermions are now possible mediators, unlike for 
d — 6. We illustrate the operator decomposition for operator #2 from Table 3, showing the 
corresponding Feynman diagram in Fig. 3. 
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Notice that only the minimal mediator content necessary to obtain each possible d = 8 
operator is shown in Tables 3 and 4. In other words, although there is always a particular 
set of exotic particles whose exchange induces at tree-level the d = 8 operators considered, 
this set might not be unique. Nevertheless, for each operator, the particle content shown in 
the tables is contained in all other possible sets of mediators leading to it. 

From both Tables 3 and 4, and from Table 2, one can easily read off the following key results 
for the operators considered : 

• There is no way to write down a d = 8 operator without involving a mediator (pin- 
pointed in boldface) which also generates d = 6 four-lepton interactions. 

• In order to build Onsi an d to cancel the dangerous (or all) NSI d = 6 contributions, 
at least two new fields are needed. 

This implies that fine-tuning - or hopefully symmetries - will be required if all d — 6 NSI 
are to be cancelled, Eq. (45) 8 . For model building, one may use the tables as follows: in 
order to generate a pure 0Nsi-type operator, it is necessary to choose effective operators 
such that Eq. (34) is fulfilled, i.e., interactions with four charged leptons are suppressed, 
and that Eq. (45) is satisfied, i.e., the NSI contributions from d = 6 operators cancel. The 
two simplest methods to build a model leading to a pure Onsi interaction are: 

1. To select from the tables those effective operators marked as Onsi- 

2. To linearly combine the effective operators in either of the tables to obtain an Onsi 
structure. One possibility is to choose any combination of at least two non-0NSi 
operators which are linearly independent in the BR basis (not considering Of 3 ^). 9 

The necessary mediators can then be directly read off from the tables; as the next step, the 
d = 6 cancellation conditions should be translated into relations among the couplings. 

Note that, in addition, there might be flavor dependent conditions and other constraints, 
which means that our tables can only serve as hints on how to build the simplest models. 
For example, one may have to worry about electroweak precision data, flavour changing 
neutral currents, non-unitarity of the PMNS matrix, loop constraints, and chiral anomalies 

8 Recall that this condition ensures that, in addition to avoiding lepton flavour violation among four 
charged fermions, other putatively dangerous d = 6 couplings are suppressed, such as for instance possible 
contributions to the very precise measurement of Gf determined from muon decay. Note as well that, in 
principle, one could avoid to impose such a strong cancellation condition by assuming very large couplings 
among the new heavy fields, and very small values for the couplings between those heavy fields and the SM 
fields which induce d — 6 operators. However, since the product between these two types of couplings will 
be present in the d = 8 operator (as in our toy model), the d = 8 couplings would be effectively suppressed 
as well and extreme fine-tuning would be needed. 

9 In short, the linear combination of two vectors involves only one free parameter (aside from the nor- 
malization). The condition in Eq. (34) amounts then to a linear equation with only one parameter, which 
can always be solved for. Since the vectors are linearly independent, they cannot cancel each other, which 
means that there will be non-vanishing NSI. If, on the other hand, one chooses linearly-dependent vectors, 
there will be no d = 8 operator at all - neither OnsIj nor the harmful one. 
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Figure 3: Example for a fully decomposed operator. The diagram corresponds to #2 of Table 3. 

if exotic fermions are introduced 10 . Also, vectorial scalar 577(2) doublets call for a deeper 
theory when present, as discussed earlier. 

Other such constraints can result from interactions of the EEEE-type, which we show 
in Table 5. Although these interactions do not produce NSI, care is mandatory when one 
introduces mediators which could induce such interactions. For example, operator #36 not 
only produces NSI, but will also lead to potential non-unitarity (through the operator in 
Eq.(lO)) and other unwanted d = 6 effects, and the operator #61 from Table 5 potentially 
leads to charged lepton flavor violation. 

For the operators in Table 3, our toy model is seen to be the only possibility using only two 
new fields, namely 2 s +1 / 2 an d 2^_ 3 / 2 . It combines operators #7, #8, #13, and #14, which 
correspond to the four effective interactions in Eq. (52) in our toy model. The table also 
allows to conclude that it is as well the most general version of the model with only two 
fields, while a simpler version might, for instance, only include #7 and #14. Recall that 
source and detection NSI cannot be created from our toy model, while matter NSI for all 
flavours are allowed. All Onsi operators obtained in Table 3 correspond to the combination 
of operators of the BR basis in Eq. (22) which are thus equivalent to Eq. (4) 

In Table 4, the simplest possibility to build a pure (9nsi and no d = 6 interaction requires 
at least three fields, namely 3q, and 3f_ 1: which may come from a large number of 
possible operator combinations. For example, one may combine operators #35 and #48. 
As discussed in Sec. 2, such a model could have correlations between source and matter NSI. 
Note that neither these models, nor our toy model, involve fermions 11 , which means that 
they cannot generate corrections to the unitarity of the PMNS matrix (through contributions 
to the operators in Eqs. (10) and (11)) nor to electroweak data (through contributions to 
the operators in Eqs. (11) and (12)), or at least not at leading order. 

10 This concerns for instance several examples in Table 3. In general, in order to cancel the chiral anomaly 
new vector-like fermions may be introduced. In the tables we just show the smallest number of mediators 
which can induce the d = 8 operators. 

11 More precisely, they do not involve Yukawa couplings linking the exotic and standard fermions. 
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(LE) (tit 2 ) (EH* ) (H t \t 2 tL) 


3/4 


1/4 




9 s , qi/-R , 9 £/-R 
^+1/2 °-l Z -3/2 


27 


(Ut 2 H*)(Y)(H t \t 2 L)(E 1p E) 


-1/2 


1/2 




11) 1 

A + l -l 


28 


(Lt[t 2 H* ) (Y)(H t \t 2 tL) (E 7p E) 


-3/2 


-1/2 




in i o L l R 


29 


(LY E c )(\t 2 ) (WH* ) ( lp )(H T \r 2 L) 


1/2 


-1/2 




Of 1 -\ L / R 1 r\L/R 

Z -3/2 X -l Z +l/2 


30 


(LYE c )(tit 2 )(WH* ) ( 1p )(H t \t 2 tL) 


3/2 


1/2 




J_ L / R 1 r\L/R 

Z -3/2 J -l Z +l/2 



Table 3: Complete list of LLEE-typc d = 8 interactions which involve two SM fields at any possible 
vertex of interaction (field bilinears within brackets). The columns show an ordinal for each operator, the 
d = 8 interaction, the corresponding combination of interactions in the BR basis, whether Onsi is satisfied 
and the necessary mediators, respectively. Those mediators leading as well to d = 6 operators in Table 
2 are in boldface. The superscript L/R indicates massive vector fcrmions. The flavor structure is to be 
understood as TP L a E 5 E^. 
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# 



Dim. eight operator 



Combination (L? L a )(L 6 H) 



a>iii 

^LLH 



A>331 
U LLH 



ri33 

^LLH 



A>313 

^LLH 



->333 sn 
<LLH U NSI 



.' r ? 



Mediators 



31 
32 
33 
34 
35 



(LYL)(L lp L)(H^H) 
(Lj p tL)(Lj p tL)(H^ H) 
(LYL)(L 1p tL)(H^tH) 
(LYtL)(L 1p L)(H^tH) 
(-ie abc )(LYr a L)x 
{L lp T b L){H^T c H) 



3V 
n 



1 1 



3u 




3i? 




Combination (L 13 L a )(L s H)(W L y ) 



3G 
37 
38 
39 
40 



(LYL)(Lh)( 1p )(Wl) 

(LYL)(LtH)(~/ p )(WtL) 
(LYtL)(LtH)( 1p )(H^L) 
(LYrL)(LH)( lp )(H^fL) 



1/2 
3/2 



1/2 
-1/2 



( 



c abc 



)(LYr a L)x 



1/2 
1/2 
1 



1/2 
1/2 
-1 



1/2 
-1/2 



/ 

/ 
/ 



(Lr b H)( lp )(Wr c L) 



Combination (L^L a )(L d H*)(L y H) 

{LYL){Ur 2 H*){ lp ){H T ir 2 L) 
{LYL)(Ltw 2 H*){ 1p ){H t w 2 tL) 
(LYtL) (Ltit 2 H* ) ( 7p ) (H t it 2 L) 
{LYTL){LiT 2 H*){ lp ){H T iT 2 fL) 



1 1 



1 n 

-,L/R 



1< 

a 



1 H 
1 H 
1 B 



oL/R 

J o 

oL/R 
J 



41 

42 
43 
44 
45 



-1/2 
-3/2 



1/2 
-1/2 



( 



If 



a be 



)(LYr a L)x 



■1/2 
-1/2 
-1 



1/2 
1/2 
-1 



1/2 
-1/2 



(LT b iT 2 77*)( 7p )(i7 i ir 2 r c L) 



1 1 



1 -C-/-R 
L -l 



1 1 



or 
°0 



+ 3^ R 



or 
a 



-1 ^ u -l 

,1,/R 
'-1 



Combination (Z^(L c ) 5 )((Z^) a i 7 

(LiT 2 L c )(I7iT 2 L)(H i 'H) 

(LTlT 2 L C )(I7iT 2 TL)(WH) 

{LiT 2 L c )(p\T 2 fL){H^TH) 
{Lt\t 2 L c )(ISw 2 L){WtH) 



4G 
47 
48 
49 
50 



)(fftjf) 

1/4 
-3/4 



-1/4 
-1/4 



(-ie a6c )(Lr a ir 2 i c )x 

(I?iT 2 T b L){H^T C H) 



1/4 -1/4 
-1/4 1/4 
-1/2 -1/2 



-1/4 
-1/4 



/ 

/ 
/ 



r 

3! 



i 
l 

3! 
3! 



Combination (h? H^)((L C ) S H)((L c ) a L-y) 



51 (Ut 2 H*){H t L c ){L c \t 2 L) 1/8 


-1/8 


1/8 


-1/8 


1/8 


/ 


l s -i 


+ 1§ - 




52 (LTiT 2 H*)(H T L c T)(I7iT 2 L) 


-3/8 


3/8 


1/8 


-1/8 


1/8 


/ 


4 


ql/S 


4- 


53 (Lfir 2 H*)(H T L c )(L^it 2 tL) 


-3/8 


-1/8 


-3/8 


-1/8 


1/8 


/ 


3 s _i 


+ lo " 


1- 3 L/ * 


54 (Ut 2 H*)(H t tL c )(IS\t 2 tL) 


3/8 


1/8 


-1/8 


-3/8 


-1/8 




3ii 4 


ql/S 
" °0 


4- 1%* 


55 (-ie a6c )(LT a iT 2 iT)x 


3/4 


1/4 


-1/4 


1/4 


1/4 




Six 4 


nL/R 
' °0 


4- i L _T 


{H T T b L c )(l?iT 2 T C L) 




















Combination (Z^(L c ) 5 )(i?t (17% 

56 {Lit 2 L c ){17H*)(H t \t 2 L) 




















1/8 


-1/8 


-1/8 


1/8 


1/8 


/ 


l'-i 


+ lo - 


\- itT 


57 (LTiT 2 L c )(17fH*)(H T iT 2 L) 


3/8 


1/8 


-3/8 


-1/8 


-1/8 




3ii 4 


oL/R 
' °0 


+ i L T 


58 {L\t 2 L c ){JJtH*){H t \t 2 tV) 


-3/8 


3/8 


-1/8 


1/8 


1/8 


/ 


l-i 4 


ni/R 

" ^0 


+ 3 L/ * 


59 (LTiT 2 L c )(I7H*)(H T iT 2 fL) 


-3/8 


-1/8 


-1/8 


-3/8 


1/8 


/ 


3ii 




h 3% R 


60 (-ie afcc )(Zr a ir 2 i c )x 


3/4 


1/4 


1/4 


-1/4 


1/4 




3ii 4 




4- 3 L/ * 


(T7 T b H*){H t it 2 t c L) 





















Table 4: Same as Table 3, but for the LLLL-type operators. Note that in this case the relationship 
between flavor structure and symbol is not unique. We show the flavor structure for each group separately. 
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# 


Dim. eight operator 


Ceeh 


Mediators 


61 




1 


1 v 


62 


(EYE)(EH T )( lp )(H*E) 


1 


1 V i nL/R 


63 


(EYE)(EW)( lp )(HE) 


1 


1 u , nL/R 


64 
65 
66 


(E 13 E c5 )(E E a E 1 )(H ii H) 
(EH^)(E C H)(WE) 
(EE C )(H^W)(EH) 


1/2 
1/2 
1/2 


IS 1 nL/R - nL/R 
L -2 -T ^_ 3 / 2 T ^+1/2 

L -2 "+" Z -3/2 Z +l/2 



Table 5: Effective d = 8 operators of the EEEE-type. The columns shows an ordinal for each operator, 
the corresponding coefficient in the BR basis and the tree level mediators, respectively. The flavour structure 
is given in the first and fourth rows. Although these operators are not NSI operators, they share with the 
latter some common mediators which may induce charged lepton flavor violation. 



E s H H E 1 H E s E 7 H H Lp L a H 




Figure 4: Examples for each category of diagrams which lead to d = 8 operators and require couplings 
of the new fields both to SM bilinears and to only one SM field. 



4.2.2 New interactions involving only one SM field 

Beyond the operators in the tables above, a much larger number of effective operators is 
obtained if, in addition to the interactions with SM bilinears, couplings between one SM 
field and two exotic fields are allowed in the fundamental theory [23]. The resulting d = 8 
operators are diagrammatically illustrated in Fig. 4 and fall in three categories, which 
contain the following SM bilinears at the external vertices: 

1. (LE)- or (LL)-type interactions with new fields. At least one of the mediators will 
necessarily induce some of the d = 6 interactions among four leptons discussed earlier 
(corresponding to the external vertices in the figure), and the couplings will thus be 
subject to the corresponding constraints. The fundamental interactions describing 
the internal vertices, however, may not be related to the previously discussed d = 6 
interactions. 

2. (LH)-tjpe interactions. In this case, the mediators do not necessarily induce any 
dangerous d = 6 operator involving four leptons, even if there are some common 
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mediators. The connections previously studied linking d = 6 and d = 8 operators 
do not need to hold. Nevertheless, these type of interactions involve exotic fermions 
(577 (2) singlets or triplets) and are constrained by non-unitary contributions to the 
PMNS matrix and some of them also by electroweak precision data, see, e.g., Ref. [23]: 
Fig. 4, center, illustrates that this class of diagrams is connected to one of the d = 6 
operators in Eqs.(lO) and (11), or a combination of them. 

3. (EH)-type interactions. These type of interactions are suggestive. The mediators 
may not induce dangerous four-fermion d = 6 operators. Furthermore, they do not 
introduce corrections to the PMNS matrix at leading order. They involve exotic 
leptons, however, which are typically strongly constrained by electroweak precision 
tests [15]. Fig. 4, right, illustrates that this class of diagrams is connected to the 
d — 6 operators in Eq. (12). 

Possible "mixed" diagrams, that is, diagrams involving two different SM bilinear couplings, 
will combine the corresponding properties. For instance, a model containing both (LE) and 
(LH) couplings to exotic mediators will simultaneously induce some of the d = 6 operators 
in Table 2 and some of the operators in Eqs. (10)-(12) which induce non-unit arity. 

It is easy to show that the vertex involving just one SM field (L, E or H) requires that the two 
exotic particles attached to it have different SU (2) x U(l) charges. Indeed, we have explicitly 
checked that all of these possibilities require at least two new fields to be phenomenologically 
viable, i.e., are not simpler than the cases discussed prior to this subsection. 

The scenarios in diagram #2 and specially #3 in Fig. 4 are appealing alternatives, as 
none of them is correlated to harmful d = 6 interactions (i.e., four charged-fermion lepton 
couplings), and #3 does not induce non-unitarity either. Furthermore, these two examples 
are Cnsi operators. Indeed, the exchange of a singlet fermion 1^ and a charged scalar 11 1 
shown in j^l gives schematically 



Here the projection onto the BR basis shows that it complies with the d — 8 cancellation 
conditions, Eq. (34). The mediator l s _ l could induce in addition d = 6 effective interactions 
if it would also couple to SM lepton doublets, as shown in Table 2, but such couplings are 
not mandatory. In contrast, the PMNS unitarity constraints should be relevant, as a singlet 
exotic fermion is involved. 

Turning now to type #3 and the scenario with an exotic doublet fermion 2^ 2 and a charged 
scalar lf_ l5 the resulting effective operator for this example is of the form 



and is thus again of the Onsi type- Furthermore, in this case the interactions neither lead 
to non-unitarity, nor any d = 6 operator in Table 2 needs to be induced if the charged scalar 
does not couple to SM lepton doublets (in other words, the d = 6 complete cancellation 
conditions in Eq. (45) can be implemented as well). Other scenarios of the kind just discussed 



(LH)(E)(E)(^L) = -\{o\ EH ) - \{ol EH ) . 



(58) 



(EH)(Z)(L)(rfE) = -\{0\ EH ) - \{Ol EH ) , 



(59) 



23 



do not necessarily have to lead by themselves to Cnsi structures: for them, cancellations 
similar to those in our toy model could be considered. However, it remains to be explored 
how difficult is to circumvent the constraints which electroweak precision tests impose on 
exotic leptons, and whether the necessary cancellations are feasible without running into 
extreme fine-tunings, for instance enlarging the scalar sector of the theory. 

During the completion of this work, Ref. [23] appeared. It explores (but is not limited 
to) the possible exchange of exotic fields which in our notation have quantum numbers 
of a scalar 11 1 (to obtain d = 6 NSI) and of a fermion 1^ (to obtain d — 8 NSI). The 
latter induces also d = 6 interactions, which lead to non-unitary contributions to the PMNS 
matrix, as it is well known and is further explored in that reference. Ref. [23] performs 
a systematic topological scan of the d = 8 operators, based on Feynman diagrams, trying 
to obtain the interaction Onsi directly from just one Feynman diagram while avoiding any 
harmful d = 6 and d — 8 contribution. Our tables correspond to the topologies 2 and 3 
in this reference, whereas the previous paragraph in this subsection would correspond to 
their topology 1. Since all possibilities in our tables contain at least one mediator leading 
to harmful d = 6 effects if one does not allow for cancellations, Ref. [23] effectively exclude 
topologies 2 and 3 in their scan (apart from our #46, which does not induce harmful d = 6 
four charged lepton interactions, but the mediator 11 1 is constrained otherwise, as we and 
Ref. [23] discussed before). Therefore, our work is complementary to that reference. Note 
that they find that the NSI in matter and the NSI at source or detector are correlated in all 
of their examples by the non-unitary effects of the heavy fermions, whereas it is easy to see 
that uncorrelated scenarios are achievable when one allows to combine different operators 
from our Table 3 (such as #7, #8, #13 and #14). As the most important difference, we 
relate the operators obtained from mediator exchanges to a complete basis of independent 
operators, which allows us to deduce the general cancellation conditions. 

5 Summary and conclusion 

In this study, we have discussed the possibility of large non-standard interactions (NSI) in 
the neutrino sector. Since any model of new physics has to recover the Standard Model at 
low energies, we have required gauge invariance under the SM gauge group and studied the 
possible effective theories. The focus is set on purely leptonic NSI, that is, on operators in 
which the only fermion fields appearing are leptons. Our analysis has been based on the 
full (analytical) decomposition of all possible dimension six and eight effective operators, 
which can be induced at tree-level by any hypothetical beyond the SM theory. Special 
focus has been set in the scenario in which the exotic mediators couple to SM bilinear field 
combinations. 

The aim is to gauge the theoretical price of achieving phenomenologically viable large neu- 
trino NSI, and to establish the minimal constraints that models have to respect for this 
purpose. Our main requirements were: 

• Interactions with four charged leptons have to be absent or highly suppressed, since 
these would lead to charged lepton flavor violation or corrections to Gf- 
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• When analyzing NSI from d = 8 operators, any d = 6 contribution among four leptons 
is not allowed or has to be very suppressed, since this would either be the dominating 
NSI (if harmless), or lead to unacceptably strong interactions among four charged 
leptons (if harmful). 

The NSI operators obtained have been expanded in a complete basis of independent opera- 
tors, which has allowed us to consistently consider cancellations among the contributions of 
different operators. This new approach has established the general cancellation conditions 
which the model parameters have to fulfill, to avoid four charged lepton interactions when 
the exotic mediators couple to SM bilinears. 

We have then studied the required complexity of any realistic model, such as what is the 
number of mediators and/or the type of cancellations needed. In short, we have demon- 
strated that is not possible to create such NSI at d — 6 for all flavour channels. For d — 8, 
we have constructed a simple toy model in which the necessary cancellations occur. It in- 
troduces two new SU(2) doublet mediators, a Lorentz scalar and a vector, which induce 
the desired large NSI without any dangerous d = 6 flavour-changing transitions among four 
charged leptons or non-unitarity corrections. It also allows for matter NSI uncorrelated 
with source or detector effects. Furthermore, we have obtained and analyzed the general 
classification of d = 8 interactions in a systematic way. 

More precisely, for the d = 6 operators it is shown at the effective operator level that, if the 
four charged lepton contributions have to be exactly cancelled, it is not possible to obtain 
e™ under the above assumptions. In addition, there are certain connections between the 
source NSI at a neutrino factory, and the matter NSI, such as = — (e* )*. We have 
shown that there is only one viable possibility for a d = 6 interaction with one mediator 
only, which is the well-known antisymmetric operator from Ref. [26]. Nevertheless, there 
are other options with more than one mediator in order to cancel all related interactions 
involving four charged leptons. Altogether, if the full decomposition of the operators is 
taken into account together with the current bounds, only might be viable with only 
one mediator and be as large as order unity. In this case, the current bounds do not require 
that the four charged lepton contributions cancel exactly. 

For the d = 8 operators, we have shown that at least two new fields are required to avoid 
the undesired d = 6 and d = 8 interactions involving four charged leptons. In fact, when the 
mediators of a d = 8 effective operator couple only to SM bilinears, there will always be at 
least one field leading as well to d = 6 contributions. These have to be cancelled in each case 
by fine-tuning or symmetries. This result can be directly seen from Table 3 and Table 4, 
which list the possible dimension eight operators, including the mediators, for fundamental 
interactions involving two SM fields. In those tables, the mediators highlighted in boldface 
lead to dimension six interactions as well. 

Our simple toy model, with two bosonic SU(2) doublets added to the SM content, leads 
directly to the desired d = 8 NSI and fulfills our minimal requirements. Notice, however, 
that we have not considered constraints from loop effects, neither for this toy model nor 
for the general operator analysis, as it is beyond the scope of the present work. These 
corrections could be specially relevant for some models when considering d = 8 interactions, 
as the one-loop corrections could spoil the d = 6 cancellation conditions. They should be 
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taken into account in any model which aims to be realistic. 

Similar considerations apply to the scenarios in which exotic couplings to both one SM field 
and SM bilinears are simultaneously allowed. In order to induce then large d = 8 NSI 
and no d = 6 couplings among four leptons, a minimum of two exotic mediators is once 
again needed. Some simple candidate models may not even require (strong) cancellation 
conditions and deserve further exploration. This is the case, for instance, when a singlet 
scalar and a fermionic doublet are added to the SM content. 

As far as the connection between source and matter NSI is concerned, we have demonstrated 
that it depends on the operators used. For example, several of the d — 8 operators in Table 3, 
or combinations of them, will only induce matter NSI, while those requiring singlet or triplet 
fermionic mediators may induce correlations (through non-unitary corrections to the PMNS 
matrix). On the other hand, all d = 8 operators in Table 4 will, in principle, allow for a 
connection between source and matter NSI independent of the mediators used. Therefore, 
it might be very well possible to detect matter NSI without source or detector effects, in 
the absence of fermions as exotic particles, such as illustrated by our toy model. Note as 
well that the models based on Table 4 require at least three new fields, which means that 
a source and matter NSI connection might be more easily achieved through non-unitary 
corrections to the PMNS matrix. 

In conclusion, we have demonstrated that the minimum complexity of a realistic model 
leading to large NSI and no charged lepton flavor violation requires at least two new fields 
inducing d = 8 NSI couplings. We have determined the possible SM charges of those medi- 
ators and the cancellation conditions for the dimension six interactions among four leptons 
that they simultaneously induce in most cases. These cancellation conditions translate into 
precise relations among model parameters. One exception might be e™ , which might be 
created at the dimension six level. Our results imply a number of constraints such that the 
observational prospects do not seem bright, specially as we did not identify some symme- 
try which would account for them. On the other side, we showed that large NSI are not 
excluded, and we found out which conditions are necessary to satisfy for any model to be 
viable. We agree that those conditions should be justified by symmetries or other arguments 
for the model to be credible. Until such justification is maybe found in some model, we 
leave it up to the reader to decide on the perspective for large NSI. 
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A On non-standard four neutrino interactions 

Although interactions among four neutrinos hardly contribute to laboratory processes, there 
has been some interest in the literature in the context of flavor oscillations in astrophysical 
environments, such as dense neutrino gases; see e.g. Ref. [91] and references therein. The 
direct laboratory bounds on these interactions are naturally extremely weak, see Refs. [92, 
93]. In this appendix, we discuss these four neutrino interactions in our gauge invariant 
framework. 



Effective operator formalism 



Since the four neutrino interactions require interactions with four lepton doublets, they only 
appear for the LLLL operators. In this case, Eq. (24) reads, including the four neutrino 
interactions, 



/ — - \ Ot^f _ 

5J? cS = — (CS^J (*VPi"«) (*YPjA) 

i/2 \ a 7 

+ ^ (ch +cl+^ 2 (ci? H + cif H - c\f H - cm ) j (pyplQ (*VPiA) 

i / 2 \ a 7 

+ h.c. . (60) 



The first point one notices is that the four charged lepton and four neutrino interactions 
share for d = 6 the same coefficient C\ L + C\ L . This means that for d — 6, any bound 
from charged lepton flavor violation etc. can be directly translated into the four neutrino 
interactions. This is illustrated here with one example. For f3 = \i and a = 7 = 5 = e, the 
bound from \i — > eee can, apart from some SU(2) symmetry breaking effects, be directly 
transferred to the four neutrino interactions. In our notation, one has 

B^e)^ (^) 2 = g, (61) 

where the non-standard parameter is defined as F = {C\ L + C| L )/A 2 - as often done 
in the literature. The current bound Br(/i — > 3e) < 10 _ 12 (90% CL) [94] then directly 
translates into F < 10~ 6 GV, which is far below any laboratory bound or even the bound 
from primordial nucleosynthesis. Of course, it is dependent on the participating flavors and 
somewhat looser for combinations involving the r, but this procedure illustrates the generic 
argument. Note that the bound for a vector mediated interaction, such as often discussed 
in the literature, turns out to be the same in this case. 

As discussed in Sec. 2, Eq. (33) should be satisfied for any realistic model in order to avoid 
these bounds. As we can read off from Eq. (60), however, the d = 6 coefficients for the four 
charged lepton and four neutrino interactions are exactly the same, which means that there 
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will not be any four neutrino interactions in that case. As a consequence, one has to go to 
d = 8 with the interactions being suppressed by A 4 . 

For d — 8, the corresponding Eq. (34) to suppress the harmful interactions among four 
charged fermions can be implemented in qualitatively different ways. For example, if Cj}} H = 
—Cf^j and C]^ H = there will be no four neutrino interactions but NSI, whereas 

for Cjjjj + Cfjjj = C]^ H + Cllh 7^ 0, there will be both four neutrino interactions and NSI. 
As it is demonstrated below, both possibilities can be realized within the model framework 
in this study. 

Model analysis 

In order to find models for large four neutrino interactions at d = 8, the same argumentation 
as in Sec. 4 is needed. First of all, Eq. (34) has to be satisfied to suppress the four charged 
lepton processes. Second, the d = 6 contributions to the NSI have to be cancelled, since 
there are strong bounds, i.e., Eq. (45) has to be satisfied. As an additional condition, one 
can not have (c/., Eq. (60)) 

CH'h + Clf H + C\f H + = (62) 

because such an operator will not contribute to the four neutrino interactions. The relevant 
decomposed operators can be found in Table 4, where one can easily read off if Eq. (62) is 
satisfied. Furthermore, note that operators which only induce C^ H will not be useful for 
the four neutrino interactions. We find from the table that operators #35, #40, #41, #43, 
#44, #46, #48, #49, #51, #52, #54, #56, #57, and #58 do not contribute to the four 
neutrino interactions. This implies that the possibility pointed out in the main text, i.e., 
to combine #35 and #48, does not lead to four neutrino interactions. One has to use more 
complicated combinations by the combination of different operators, such as #32 and #50 
to satisfy Eq. (34), and #48 (which satisfies Eq. (34)) to introduce an additional mediator 
to cancel the d = 6 NSI. Then the four neutrino interactions can be constructed with three 
different mediators, where only #32 and #50 contribute to the four neutrino interactions. 
Constructions with less mediators are, under the assumptions in this study, not possible, 
which is different from the NSI, which can be generated from two mediators. 

As soon as a specific model is known, the relationship among source and production NSI, 
matter NSI, and four neutrino interactions can be easily calculated using Sec. 2 and Eq. (60). 

In summary, for the d = 6 four neutrino interactions, gauge invariance implies that they face 
the stringent bounds from charged lepton flavor violation, such as from /i to three electrons. 
Therefore, large four neutrino interactions have to come from d = 8 effective operators. 
From the model point of view, having four neutrino interactions is even more complicated 
than having large NSI, since at least three different mediators are needed in the framework 
discussed in this study. 
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